Wounded quarks and diquarks 
in high energy colhsions 



A. Bialas and A. Bzdak* 
M. Smoluchowski Institute of Physics 
Jagellonian University, Cracow^ 

February 26, 2008 



Abstract 

Particle production in Au — Au, Cu — Cu, d — Au andp— p collisions at 200 
GeV cm. energy are analyzed in the wounded quark-diquark model. Existing 
data are well reproduced. Emission functions of wounded and unwounded 
constituents are determined. Implications for the collective evolution of the 
system are discussed. 



1 Introduction 

It is now well established that the idea of a "wounded" source of particles [1] turned 
out rather useful in description of particle production from nuclear targets [2]. It 
seems therefore interesting to verify to what extent the recent high energy data 
confirm this hypothesis and what are its consequences for the evolution of the system 
created in heavy ion collisions at high energy. This is the subject of the present 
paper. 

A wounded source, by definition, emits a fixed density of particles, independently 
of the number of collisions it underwent inside the nucleus. To understand the 
physical meaning of this concept, let us recall the original argument leading to this 
idea. 

The idea can be understood if one observes that the process of particle production 
is not instantaneous. This was first noted (and the consequences for scattering at 
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high energies derived) by Landau and Pomeranchuk [3] . A simphfied version of the 
argument can be presented as follows. 

Consider a particle created in a high-energy collision. In the reference frame 
where the longitudinal momentum of this particle vanishes, the minimal time nec- 
essary for its creation is to ^ 1/""^* where rrit = \/ m? + pi is its energy. Consider 
now the "laboratory" frame where the target nucleus is at rest. In this frame the 
particle in question acquires some longitudinal momentum, the time is multiplied 
by the Lorentz factor, and we have 

. . E coshyiab 

^ = 7^o> — = , 1) 

mf rrit 

where E is the energy of the particle. Consequently, the uncertainty of the distance 
from the collision point to that at which the particle is created (i.e. the resolving 
power of the longitudinal distance) is 

l = vt>'^^. (2) 
rrit 

When this distance is longer than the size of the nucleus, Z{b), at a given impact 
parameter i.e. when the rapidity of the produced particle is large enough, the 
particle cannot resolve individual collisions and therefore it is natural to admit that 
its creation may be insensitive to the number of collisions of the source. This is 
the origin of the idea of wounded sources and it is clear that it makes sense only 
when one considers production of particles with the laboratory rapidity exceeding 
that determined by the condition / > Z{b). But at small laboratory rapidities 
particle production from a fast-moving source is anyway expected to be small (if 
not negligible). Therefore the concept of the wounded source can be applied in 
practice in the whole rapidity region. 

It was recently shown [4] that the idea of wounded constituents of the nucleon 
can successfully explain the centrality dependence of particle production at 77 = in 
Au — Au coUisions at RHIC [5] . Assuming that high energy interactions of nucleons 
are dominated by independent interactions of its two constituents, a quark and a 
diquark (each one producing the same particle density), it was possible to explain 
the small momentum transfer elastic pp and vrp scattering [6] as well as energy and 
centrality dependence of particle production in Au — Au collisions at the same time^. 

Encouraged by these results, in the present paper we extend the analysis to 
rapidities in a wide range beyond the central region using data on particle production 
in p — p, d — Au, Au — Au and Cu — Cu collisions at 200 GeV cm. energy. 

Our main conclusion is that the data are in reasonable agreement with the 
wounded quark-diquark model in a wide range of rapidities excluding, however, the 
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narrow target and projectile fragmentation regions. In the fragmentation regions 
two additional contributions could be identified: 

(i) First, it is necessary to account for the decay products of the constituents 
which did not interact (" unwounded" ) but belong to a wounded nucleon. Since such 
constituents are coloured, they have to decay into observable particles. 

(ii) One also has to take into account particle production by secondary interac- 
tions of particles produced by the projectile inside the nucleus. Since, as already 
explained, the fast particles (in the rest frame of the nucleus) are produced outside 
the nucleus, the secondary interaction may only produce particles which are slow in 
the rest frame of the nucleus. 

In the next section we present a general formulation of the model. The d — Au 
collisions are discussed in Section 3 and the nucleus-nucleus collisions in Section 4. 
Our conclusions are listed in the last section where also some comments are included. 



2 General formulation of the model 

In the central rapidity region the model predicts that the density of particles pro- 
duced in collision of nucleus A with nucleus B is given by [4] 

^^^^ = PAsiy) = w^^F^iy) + w^^Fj,{y), (3) 

where w^^^ is the number of wounded constituents (quarks and diquarks) in the 
nucleus A{B), whereas Fa and Fb are particle densities emitted by one wounded 
quark or diquark in the nucleus A and S, respectively. Evaluation of w^'^^ for various 
processes considered in this paper is given in the Appendix. 

In the cm. system we have, of course 

FB{y)^FA{-y) = F{y), (4) 

where from now on, y shall always denote rapidity in the cm. frame. 

As explained in introduction, to extend the description to the region close to the 
maximal and minimal rapidities, it is necessary to take into account at least two 
other contributions to particle production. Taking this into account, the Eq. (3) is 
modified into 

PAB{y)=w^;^^F{y)+w'~^^F{-y) 

+ wfu{y) + wfu{-y) + Cs(y) + C^(-|/), (5) 

where U{y) = Usiy) = UA{—y) represents the contribution from the decay of one 
unwounded constituent (but belonging to a wounded nucleon) and w^ab num- 
ber of unwounded constituents in the nucleus A{B). It can be simply determined 
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using our assumption that each nucleon consists of two active constituents (a quark 
and a diquark) 

Wa,B + '^A,B = '^^A,B, (6) 

where wa,b is the number of wounded nucleons in the nucleus A{B). 

The terms Ca and Cb represent the "intranuclear cascade", i.e., the contribution 
from secondary interactions of particles created inside the nucleus. 

The contributions Ua and Ca are expected to be significant only at y close to 
—Y (the cm. rapidity of the nucleon in nucleus ^4). Similarly the Ub and Cb are 
expected to be significant only at y close to Y. 

In nucleon-nucleon collisions, the " cascade" contribution is absent and we obtain 

PNNiy) = [Fiy) + Fi-y)] + w^^\uiy) + C/(-y)], (7) 

where from the analysis of elastic pp data we have found [4] that = 1.187. 
Following (6) we obtain w^^'^ = 0.813. 

The two equations (5) and (7) describe the relation between particle production 
in nucleon-nucleon, nuclcon-nuclcus and nucleus-nucleus collisions implied by the 
model. They describe the data in terms of three unknown functions. Since these 
functions have a well defined physical meaning, the important result of our inves- 
tigation is not only to verify the validity of (5) and (7) but also determination of 
the emission function of a wounded constituent F{y), the decay distribution of an 
unwounded constituent U{y) and the "cascade" contribution C{y). 

3 d-Au collision 

The PHOBOS data on charged particle production in d — Au collisions at ^/s — 200 
GeV [8] cover the cm. pseudorapidity range \ri\ < 5.3. 

The first step is to check if the data are consistent with the model at 77 = 0, 
where the simple formula [implied by (3) and (7)] 

Pd-A.(0)-PiViv(0) '^^'"^J^'\ (8) 

holds with Pnn{0) being the density of particles produced in a single pp collision 
at 77 = 0. In the further calculations we take Pnn{0) = 2.3llo!i6 [9]- Fig. 1 
dNd-Au/ drj at = is plotted versus w^+Au = wau + Wd (number of wounded 
nucleons in both colliding nuclei) and compared with the model prediction (8). 
One sees that, within substantial experimental errors, the data are consistent with 



4 



the modeP. It is also seen, however, that the data for most central collisions are 
somewhat below the results given by (8). 
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Figure 1: Wounded quark-diquark model compared with the midrapidity d — Au 
data from PHOBOS coll. 

To investigate the rapidity region outside y = 0, we first ignore the "intranuclear 
cascade" contributions Cd and Cau- This simplification allows to investigate the 
data in terms of only two functions: F{y) = = y) and U{y) = Ud{y) — 

UAu{—y), being the densities of particles produced from one wounded or unwounded 
constituent, respectively. This condition is well justified for the deuteron since the 
intranuclear cascade in this case is negligible. On the other hand, we expect a 
discrepancy in the Au fragmentation region. This discrepancy provides a measure 
of the effect. 

We construct F{ri) and U{ri) using the PHOBOS d — Au data [8]. First, to avoid 
propagation of the slight difference between the model and the data observed at 
?7 = (c.f. Fig. 1), we adjust data to the model at this point, multiplying at each 
rapidity by a constant (for each centrality) factor P^a-Aui^) I p1!^-aJS^)- '^^^^ allows 
to discuss the shape of the distributions in rapidity independently of the absolute 
normalization (which is also consistent with the model, as seen in Fig. 1). 

Next we verified that in the pscudorapidity range I77I < 3.5 the contribution from 
the decay of unwounded constituents is not necessary for description of data and 
thus one can take U {rf) = 0. It means that the data in this region can be described 
solely in terms of wounded constituents. This was confirmed by analysing the data 
on the deuteron side {0 < 1] < 5.3) in terms of three functions: Fd{i]) = F{r]), 

^The errors in model prediction reflect the inaccuracy in the pp data i.e. pjvjv(O). 
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Fau{v) = Fi—ff) ^iid Ud{r]) — U{r)). Prom this analysis we have also found that 
Fau{v) = for 77 > 3.5. 

Using these findings and the 20 — 40% centrality PHOBOS data^ it was possible 
to determine F{ri) for \ri\ < 3.5. For rj > 3.5 we used the 20 — 40% and 40 — 60% 
centrality PHOBOS data to obtain both F{ri) and f/(?7). 

In Fig. 2 the functions F{r]) and U{r]) are shown. One sees that, as expected 
[10-12], U{r]) is confined to the region close to the maximal rapidity. One also 
observes that F{ri) is negligible for r] < —3.7. 
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Figure 2: The densities of particles produced from one constituent at ^/s — 200 
GeV. Squares: wounded, F{ri); Crosses: unwounded, U{r]). 

Using these F{ri) and U (77) one could evaluate the predictions of the model for all 
centralities and in the full rapidity range. The results are compared to the PHOBOS 
[8] and STAR [13] data in Fig. 3. One sees that the agreement is satisfactory for 
77 > —4. The discrepancy for 77 < —4 can be attributed to the additional "cascade" 
contribution which is expected in the Au fragmentation region. Large errors do not 
allow, unfortunately, to perform a more quantitative analysis of this phenomenon. 

All in all, the results presented in Fig. 3 show that the model describes correctly 
the data. 

Since the functions F and U are now determined, it is possible to construct the 
model prediction for pp collisions. In Fig. 4 this is confronted with preliminary 
pp PHOBOS data [9]. We see that the agreement is very good. For comparison 

^For the 77 < one observes a certain difference between the STAR and PHOBOS data for 
— 20% centrality. Therefore we preferred to take the 20 — 40% centrahty where this difference is 
much smaller. 
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Figure 3: Wounded quark-diquark model compared to d — Au data from PHO- 
BOS and STAR coll. The data are normalized to the model prediction at 77 = 0. 
One sees a generally good agreement except in the Au fragmentation region where 
the contribution from secondary interactions inside the target apparently becomes 
important. 




Figure 4: Wounded quark-diquark model compared with the preliminary pp PHO- 
BOS data. For comparison we also present the prediction without contribution from 
unwounded constituents. One sees that it is necessary to account for the data. 
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we also present the prediction of the model without contribution from unwounded 
constituents. It appears that this contribution is indeed needed. 

We feel that this result indicates that the physics of particle production in pp 
and d — Au collisions is basically the same (apart from the trivial difference in the 
form of intranuclear cascade). 

It may be interesting at this point to speculate about the possible mechanism 
which may lead to the results for F{7]) and Uij]) shown in Fig. 2. The very broad 
rapidity distribution of particles emitted by the wounded source suggests a colour- 
exchange mechanism. The simplest possibility is the multi-gluon exchange between 
the projectile and target (every gluon exchange changes the colour of the source 
and thus may lead to particle production) . Since every constituent consists of many 
partons and since the probability of exchanging a gluon between two partons does 
not depend on their rapidity difference (gluon is a particle of spin 1), one may fairly 
easily obtain a rather broad rapidity distribution of the produced particles. As dis- 
cussed in [14, 15] this picture can naturally accommodate the asymmetric shape of 
F{r}) (the shoulder seen at 77 = —1 is -most likely- the kinematic effect related to the 
transformation from rapidity to pseudorapidity) . On the other hand, particle pro- 
duction from the unwounded source can be understood as resulting from the decay 
of a string spanned between the wounded and unwounded constituent. The details 
of this picture and its derivation from fundamental theory represents an interesting 
problem which, however, goes beyond the scope of the present investigation. 

4 Au-Au and Cu-Cu collisions 

In Fig. 5 we compare the predictions of the wounded quark-diquark model with the 
PHOBOS data [5] in Au—Au collisions. This observation shows that the intranuclear 
cascade is not effective for central nucleus- nucleus collisions'^, thus suggesting that 
the secondary interactions in the nucleus lead to particle production mainly if they 
happen on the spectator nucleons. It would be interesting to investigate this point 
in more detail when more precise data are available. 

The values of the numbers of wounded constituents for various centralities are 
given in the Appendix. 

We have found in the previous section that F{ri), the contribution form wounded 
constituents, largely dominates the spectrum in the region \ri\ < 4.5. It then follows 
from (5) that, for symmetric nucleus-nucleus collisions, the ratio of particle densities 

closer look shows that there is some indication of discrepancy in the fragmentation region 
in case of most peripheral collisions. 
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Figure 5: Wounded quark-diquark model compared with the Au — Au data from 
PHOBOS coll. 

produced by various nuclei must be a constant, independent of rj. Indeed, we have 



PAuAuiv) ^ 
PCuCuiv) W^nl 



Rau/Cu{v) = = — = Rau/Cu{Q)- (9) 



This consequence of the model is known to be very well satisfied by Au — Au and 
Cu — Cu data [16]. 

In Fig. 6 -RyiM/Cu(0) evaluated from the model is shown versus the number of 
wounded nucleons (the same for Au and Cu). One sees that it is very close to 1, 
also in agreement with data [16]. 



5 Conclusions and comments 

Our conclusions can be formulated as follows. 

(i) We have compared the wounded quark-diquark model with the PHOBOS 
data on particle production in p — p, d — Au., Cu — Cu and Au — Au collisions 
at 200 GeV cm. energy in the full range of rapidity. The shape of the rapidity 
distribution is well reconstructed. The overall normalization of the d — Au data 
is not described so well, although it is consistent with the model within the (still 
rather large) experimental uncertainties. 

(ii) The particle emission function from one wounded constituent is determined. 
Its most important features are: (a) a large maximum in the forward direction; (b) 
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Figure 6: The ratio of density of particles produced at midrapidity in Au — Au 
collisions to that produced in Cu — Cu collisions at the same number of wounded 
nucleons. For comparison we also present the results at 19.6 GeV cm. 

a significant fraction of particles emitted in the backward hemisphere; (c) strong 
suppression of particle emission in the target fragmentation region. 

(iii) Particle emission from unwounded constituents contributes only to the frag- 
mentation region of the projectile. 

(iv) Particle production from secondary interactions inside the nucleus are con- 
fined to its fragmentation region, in agreement with theoretical expectations. These 
secondary interactions seem to have stronger effect on particle production in nucleon- 
nucleus than in the central nucleus-nucleus collisions. 

(v) For the good agreement of the model with data it was crucial to accept that 
each wounded constituent emits the same density of particles. This somewhat sur- 
prising result indicates that intensity of particle emission is mainly determined by 
the colour content of the source (the colour content of quark and diquark are the 
same). This suggests the colour-exchange models as a possible mechanism respon- 
sible for particle production, in harmony with the discussion at the end of Section 



Following comments are in order. 

(i) Contrary to the common prejudice, good agreement of the wounded con- 
stituent model with Au — Au and Cu — Cu data is not in contradiction with the 
collective phenomena observed in heavy ion collisions. It only shows that in all 
hadronic collisions the early stage of the particle production process can be under- 
stood as a superposition of contributions from hadronic constituents. This does 
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not preclude further collective evolution of the system which obviously must be 
more visible in the system produced in collision of two heavy nuclei than, e.g., in 
nucleon-nucleon collision. 

(ii) The wounded constituent picture has, however, important consequences for 
the properties of the hot matter created in heavy ion collisions. First of all, it 
implies that most of the entropy must be produced already at the very early stage 
of the collision. It also implies rather early equilibration in transverse direction, 
presumably already at the level of nucleon-nucleon collisions [17-19]. Finally, the 
absence of the visible effects of the longitudinal pressure in Au — Au coUisions (which 
would modify somewhat the rapidity distribution) suggests that the early evolution 
of the plasma may be dominated by the purely transverse hydrodynamic expansion 
while the longitudinal evolution is described by free-streaming. This possibility 
(which solves the "problem of early equilibration") was investigated recently [20] 
and shown to be in agreement with data on elliptic flow. One may thus consider 
our result as an indirect confirmation of the hypothesis formulated in [20]. 

(iii) The model describes particle production in nucleon-nucleon collisions as 
emission from two (left- and right-moving) sources, both populating most of the 
available phase-space. This picture suggests presence of specific long-range, forward- 
backward correlations. It would be interesting to study these correlations in more 
detail. 

(iv) In the present investigation we have entirely neglected the possible depen- 
dence on transverse mass, thus implicitly assuming that rrit distributions do not vary 
significantly with the varying centrality of the collision. This is approximately cor- 
rect for the bulk of created particles (mostly pions at small transverse momentum) 
but it would be certainly interesting to investigate limits of the wounded constituent 
model at higher transverse momenta and for heavy particles. 

(v) The model provides a natural qualitative explanation of the "stopping" of 
a high-energy nucleon in a collision with the heavy nucleus [21]. Indeed, since the 
contribution from an unwounded constituent dominates the very end of rapidity 
phase-space and since the number of unwounded constituents in the nucleon passing 
through a nucleus is smaller than that in a nucleon-nucleon collision, one expects less 
high-energy nucleons in the former case. A quantitative estimate of this effect may 
provide information about the momentum distribution of the constituents (quark 
and diquark) inside the nucleon. 
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A Appendix: Estimate of the number of wounded 
constituents 

Al. Au-Au 

Evaluation of the number of wounded quarks and diquarks in Au — Au 
collision is described in detail in [4]. For the three centralities, — 6%, 15 — 25% 
and 35 — 45%, presented in Fig. 5 wc obtain w^^l = 320, 175 and 76, respectively. 
The number of unwounded constituents is determined from (6) with the values 2wau 
taken from [5] {2wau = 344, 200 and 93, respectively). The same procedure is used 
for Cm — Cm collisions. 
A2. d-Au 

Evaluation of the number of wounded quarks and diquarks in d — Au system is 
somewhat more complicated. 

First we assume that each wounded nucleon in Au is hit only once^, so the 
number of wounded constituents in Au is given by: w^^l = wauW^n^ where wau is 
the number of wounded nucleons in Au (provided by the PHOBOS coll. [8]) and 
wJv^ = 1.187 [4]. 

To estimate the number of wounded quarks and diquarks in the deuteron let us 
first consider the nucleon-nucleus collision. 

The average number of wounded quarks w^^\ in the nucleon which underwent k 
inelastic collisions is given by the straightforward counting of probabilities 

where Pg is the probability for a quark to interact in a single pp collision. In [4] we 
have found that Pg ^ w^''/3 = 0.395. One can write analogous formula for w^*^'' i.e. 
number of wounded diquarks in the nucleon which underwent k inelastic collisions. 

Thus, the number of wounded constituents (quarks and diquarks) in the nucleon 
which underwent k inelastic collisions is 

^{;?. = 2-(l-p,)'-(l-Pd^.)', (11) 

where pdiq is the probability for a diquark to interact in a single pp collision and 
Pdiq ~ 2pg = 0.79. 

^This is a good approximation since there is no significant difference between the number of 
wounded nucleons in Au and the number of colhsions [8]. 
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For the average number < A; > of collisions per one wounded nucleon in deuteron 
we take < k >= ncou/wd, using the values of Ucou (number of collisions) and 
(number of wounded nucleons in deuteron) given by the PHOBOS coll. [8]. 

Replacing in (11) A; by < k > we obtain^ ^^''iv''<fc> finally, the number of 

wounded constituents w'^^ in deuteron for a given Hcou and Wd reads 

= ^c^«^SU>- (12) 

This procedure gives: w^"^^ = 3.95, 3.66, 3.1, 2.3, 1.6 for - 20%, 20 - 40%, 40 - 
60%, 60 - 80% and 80 - 100% centrahty, respectively. 
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